Magnetism: the Driving Force of Order in CoPt. A First-Principles Study. 
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CoPt or FePt equiatomic alloys order according to the tetragonal Llo structure which favors their 
strong magnetic anisotropy. Conversely magnetism can influence chemical ordering. We present 
here ab initio calculations of the stability of the Llo and LI2 structures of Co-Pt alloys in their 
paramagnetic and ferromagnetic states. They show that magnetism strongly reinforces the ordering 
tendencies in this system. A simple tight-binding analysis allows us to account for this behavior in 
terms of some pertinent parameters. 

PACS numbers: 71.15.Nc,71.20.Be,75.50.Cc,61.66Dk 



Magnetism and chemical ordering are frequently cou- 
pled in alloys. On one hand, a strong magnetocrystalline 
anisotropy characteristic of equiatomic binary alloys like 
FePt or CoPt is known to be due to an Llo order that 
alternates pure planes along the [001] direction [l|-|5|. On 
the other hand, magnetism can influence phase stability 
and chemical ordering [gl, 0- Typical examples include 
the stability of the bcc a phase of iron [8] and the phase 
diagram of FeCo [1, [l^ . In the case of NiFe alloys in- 
volving neighboring late transition elements, magnetism 
plays an important role in the stability of ordered phases 
in a rather notable fashion [ll|, [l2[ . It is suspected that 
this should also be true for FePt and CoPt alloys and, 
in the case of CoPt, diffuse scattering as well as nuclear 
magnetic resonance experiments have established a clear 
relationship between magnetism and short range order 

Furthermore, these alloys, in the form of nanometer- 
sized grains, are ideal candidates for high density mag- 
netic storage applications (provided that the Llo ordered 
state is preserved). The systematic modeling of their 
thermodynamic properties, of the critical ordering tem- 
perature in particular Q, necessitates the development 
of multi-scale methods involving effective interatomic po- 
tentials. These potentials should include all the relevant 
physics of the alloy, while remaining simple enough to 
allow simulation of real life situations with hundreds or 
even thousands of atoms l3-l3]- It then becomes im- 
perative to determine whether magnetism governs ener- 
getic properties of bulk alloys and nanoalloys. Should 
this prove to be true, all interatomic potentials would 
need to include a magnetic term. Clearl y, t his is not the 
case of the potentials in current use |20M24j . 

The purpose of this Letter is to quantify the role 
of magnetism in the formation energies of Co-Pt alloys 
through non-magnetic and magnetic ab initio electronic 
structure calculations; thus providing tangible grounds 
for a magnetic interatomic potential. In the bulk form, 
three ordered phases are known to exist in the Co-Pt sys- 
tem at low temperature, corresponding to the stoichio- 
metric concentrations of CosPt, CoPt, and CoPts [25[. 



The ordered phase of CoPt is of the tetragonal Llo type. 
CoaPt and CoPts are both LI2 ordered phases of cubic 
symmetry. The bulk order-disorder transition of CoPt is 
equal to 1098 K 0. 

We carried out spin-polarized calculations in the 
framework of the Density Functional Theory (DFT) us- 
ing the ABINIT code [23 with the Generalized Gradi- 
ent Approximation (GGA) exchange correlations func- 
tionals. Core and valence electrons were represented by 
a plane wave basis and the projector augmented wave 
(PAW) potentials The adopted valence electronic 

configurations for Co and Pt are Sd^As^ and 5(fi6s^ re- 
spectively. All plane waves with energies below the cut- 
off energy were included in the basis set. The cut-off 
energies (16 eV, 22 eV, 22 eV for Co, Pt and CoPt, re- 
spectively) were chosen 25% larger than the largest de- 
fault cut-off of the element-specific potentials. Integra- 
tions over the Brillouin zone are based on a 20 x 20 x 20 
Monkhorst-Pack 2D grid which is sufficiently fine to en- 
sure the numerical convergence of all the calculated prop- 
erties. All the structures were fully relaxed using the 
Broyden-Fletcher-Goldfarb-Shannon minimization. The 
cold smearing method was used for the Brillouin zone 
integration with a smearing parameter of 7xl0~^ Ha 
leading to formation energies converged to within 10~^ 
eV/at. The calculations were performed at zero pressure; 
the relaxation of the atoms and the shape of the simu- 
lation cell are considered using the conjugate gradient 
minimization scheme. The atomic positions are relaxed 
until the forces on the atoms are reduced to within 10~^ 
Ha/Bohr. 

We first verified that face centered cubic (fee), body 
centered cubic (bcc) and hexagonal close packed (hep) 
structures were correctly reproduced in elemental Co and 
Pt. The relative stability of these various phases, as well 
as the influence of magnetism on the system can be deter- 
mined from energy versus atomic volume curves plotted 
in Fig. [1^ for Pt and Fig. [TJd for Co. Concerning Co, 
the non-spin polarized calculations predict a fee ground 
state while spin polarized calculations correct this and 
reproduce the experimentally stable hep phase (see Fig. 
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FIG. 1: Total energy as a function of atomic volume for (a) Pt 
and (b) Co as predicted by GGA calculations (open markers 
: non magnetic, full markers : ferromagnetic). For Co, the 
variation of the magnetic moment is also displayed. 



[TJd). Such behaviors of non magnetic calculations have 
been highlighted in previous calculations [29j . Regarding 
the magnetic moment, as expected, it increases when the 
lattice is expanded and vanishes when it is reduced. The 
total energies of the ferromagnetically ordered bcc and 
fee structures are slightly higher than the hep ones, 0.12 
eV/at and 0.01 eV/at excess energies per atom, respec- 
tively. Since the energy difference between fee and hep 
is relatively small and since the CoPt alloys are cubic or 
tetragonal according to the phase diagram, cobalt will be 
considered to be in its fee phase in what follows. 



TABLE I: Physical properties of Co and Pt in the FCC struc- 
ture (GGA calculations). Lattice parameter (ao) and mag- 
netic moment per atom (mtot)- DFT and experimental (when 
they exist) values are presented in brackets. 



ao 



(A) 



mtot (ij^b) 



Co 

CosPt 
CoPt (Llo) 

CoPt (Lli) 

CoPt (A2B2) 

CoPts 

Pt 



3.52 (3.53 [30]) 
3.66 (3.66 [30j) 
a=3.81 (3.81 [30j) 
c/a=0.976 (0.973 [30j) 

a=3.70 (3.80 [33]) 
b/a=1.016 (1.017 [33]) 
3.71 
3.86 (3.83 [30]) 
3.97 (3.92[34]) 



1.64 (1.62 [31]) 
1.43 (1.45 [32]) 
1.14 (1.20 [6]) 

1.08 

1.12 
0.73 (0.67 [6]) 




As shown in Table HI our results for physical and mag- 
netic properties for the fee structure are in agreement 
with both previous ab initio, and when relevant, exper- 
imental data. We have calculated a magnetic moment 
within the LDA and GGA approximation equal to 1.54 
/i^/at and 1.64 /i^/at, respectively. The LDA underes- 



timation of the magnetic moment is clearly connected to 
the well-known underestimation of the lattice parameters 
generally obtained with LDA. Conversely, the usual over- 
estimation of the lattice parameter by the GGA is not 
observed in the case of Co, leading to a good agreement 
between the experimental and the GGA determination of 
the magnetic moment. This leads us to prefer the use of 
GGA to model the CoPt system, even though Pt struc- 
tural properties are slightly better reproduced by LDA. 
Regardless of the approximation, note that all structures 
for Pt are found to be non magnetic (NM), the ground- 
state structure being NM fee as shown in Fig. [T^. 

We now focus on ordered LIq and LI2 alloys which are 
known to exist in the CoPt system. The LIq phase ex- 
hibits uniaxial anisotropy with lattice parameters a and 
c about 3.81 and 3.71 A respectively. As was the case for 
elemental Co and Pt, the physical properties obtained for 
both alloyed structures are in excellent agreement with 
experimental data and other ab initio results, as shown 
in Table [D Moreover, we have also considered a simplified 
geometry where c = a; the magnetic moment was found 
to differ by 2 % and the total energy by 0.02 eV/at. This 
suggests that the precise value of the c/a ratio is irrel- 
evant when the foci are the general physical properties 
of the structure. In particular, its magnetocrystalline 
anisotropy is principally governed by the (anisotropic) 
chemical ordering [35| more than by the deviation from 
unity of this ratio. Therefore, in the following we assume 
that c = a. 
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FIG. 2: Atom-projected electronic densities of states of CoPt 
Llo: NM (a) and FM (b). 



In Fig. O we show the atom-projected electronic densi- 
ties of states (PDOS) for the Llo-non magnetic (NM) and 
Llo-ferromagnetic (FM) phases of CoPt. As expected, 
the Pt band is wider than the Co band, about 8 eV and 
5 eV respectively. The bands above — 7 eV are predom- 
inantly of d character and also show a large hybridiza- 
tion, in this case between the Co 3d and Pt 5d states. 
In the paramagnetic phase, the Fermi level falls within a 
Van Hove peak, which, according to the Stoner criterion, 
suggests a possible instability towards a ferromagnetic 
state. Such is the case, in Fig. [2l3, where a quasi-rigid 
shift between the majority (spin-up) and the minority 
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(spin down) bands is observed. This alloy is a strong 
ferromagnet with a full majority spin band as was the 
case for elementary fee or hep Co. About 1 eV above 
the Fermi level, we observe a hybridization between Co 
and Pt d states, leading to a small induced local mag- 
netism on Pt, about 0.3 /i^/Pt atom. The PDOS for 
the LI2 phase shows roughly the same behaviour. These 
calculations are consistent with previous works (36l-[38|. 

Similar studies were conducted in the alloyed phase at 
varying concentrations of Pt (Coi_a^Pta^). The ordered 
phases LI2 (CoaPt or CoPts) have been considered, as 
well as other equiatomic phases : Lli, and the A2B2 
phase fu] . Our results are presented in Table |T1 For the 
Lli structure in particular, the b/a ratio has been found 
equal to 1.016, closely comparable to the value 1.017 of 
previous ab initio calculations [33|. The total magnetic 
moment expressed per atom as in table I decreases as 
a function of the Pt concentration: this is due to the 
fact that the magnetic moment of Pt atoms is negligi- 
ble, whereas the magnetic moment of Co atoms is almost 
independent of the Pt concentration, and even increases 
slightly with it, as already noticed [32|. 

Finally the enthalpies of formation {AH, at T = OK) 
of these alloys were determined: 



AH 



^t'tt°^(^Co,mPt) -n£;i 



Co 
fee 



l(n- 



m) 

(1) 

where E^^"^^ {nCo, mPt) is the total energy of the mixed 
Co + Pt system containing n Co atoms and m Pt atoms, 
and represents the energy per atom of the elemental 
form X (X= Co or Pt) in the appropriate reference state. 
Thus, for the non spin polarized calculations, the refer- 
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FIG. 3: Enthalpy of formation at OK as predicted by GGA 
calculations : non magnetic and ferromagnetic. 

ence states for the two elements are the non magnetic 
structures whereas for the spin polarized calculations we 
consider the ferromagnetic state for Co and the non mag- 
netic state for Pt. 



Fig. [3] illustrates the importance of magnetism in as- 
sessing the stability of the ordered phases of the Co-Pt 
system. In non spin polarized calculations, AH is found 
to be very weak for any concentration. This result is 
incompatible with experimental observations of ordered 
structures having critical temperatures above 1000 K. In 
the spin polarized calculations, AH takes negative val- 
ues between -0.07 and -0.10 eV/at for the LI2 and LIq 
phases. These phases are therefore clearly stabilized by 
magnetic effects. The Lli and ^42^2 equiatomic struc- 
tures are also clearly stabilized by magnetic effects but 
are less stable than the LIq phase. The effect of mag- 
netism is more pronounced for the Lli phase than for 
the A2B2 and LIq phases, but not sufficient to modify 
the relative stability of the phases in comparison with 
the non spin polarized calculations. The energy differ- 
ences between these three phases show that effective pair 
interactions beyond first neighbors are not negligible [llj 
and that our ab initio calculations are successful in pre- 
dicting Llo to be the most stable phase. 

Our results for the Lip phase agree also with those ob- 
tained by Alam et al. [39|; these authors (who did not 
considered other equiatomic structures) found a forma- 
tion enthalpy about -0.09 eV/at for the ferromagnetic 
state, to be compared to our value about -0.10 eV/at 
and the experimental value about -0.13 eV/at [4o|. Their 
value for the paramagnetic state is slightly smaller than 
ours: -0.03 eV/at instead of -0.02 eV/at. It should 
be mentioned that the uncertainties due the approxima- 
tions made (LDA or GGA, pseudopotentials, KKR-CPA 
in the case of disordered states) are probably much larger. 
Recent calculations {41^ also agree on the order of mag- 
nitude of the formation enthalpy (about -0.1 eV/at for 
LIq ferromagnetic CoPt) but the role of magnetism is 
not discussed. 

Even though ab initio calculations are efficient in de- 
termining a quantitative description of an alloy with rel- 
atively high precision, general trends are best understood 
using simple models such as those provided by the tight- 
binding approximation [ll|, ll2|. In the simplest model. 



the main alloy parameter is the diagonal disorder param- 
eter, i.e., the difference in atomic levels Sd = — ej[ 
compared to an average bandwidth. In this context, the 
main outcome of tight-binding studies is the justification 
of an Ising Hamiltonian for order phenomena in transi- 
tion metals [42|. In addition, effective pair interaction 
between first neighbors {V) which dominates this Hamil- 
tonian can be calculated through the Coherent Potential 
Approximation (CPA) for disordered alloys |11|. 

Fig.Hla illustrates a typical variation of ^ as a function 
of the average number of (i-electrons Nd present in the 
system for an equiatomic alloy. Although for a real al- 
loy Nd and the atomic concentration are related through 
Nd = cN^ + (1 - c)7Vf for an A^^i-c alloy, where N^ 
is the number of d-electrons (X=A or B), it is convenient 
here to consider them as independent variables. Ordering 
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FIG. 4: Typical variation of V as a function of the band-filling 
Nd for an equiatomic composition. V > means an ordering 
tendency, (a) 6d— 2 eV, dnd = eV and (b) 6d— 2 eV, dnd — 
3 eV. Schematic representations of band filling for CoPt (c) 
non magnetic and (d) ferromagnetic. 



ments on Co that matters, more than the type (ferro- 
magnetic, antiferromagnetic) of long range order which 
is stabilized at low temperature. This is clear when us- 
ing the disordered local moment (DLM) picture where 
local moments interact through effective pair interactions 

To summarize, from the analysis of our ab initio cal- 
culations within a simple tight-binding model, we have 
shown that order in CoPt alloys, in particular at the 
equiatomic concentration, is principally driven by mag- 
netic effects. Off-diagonal disorder plays a role by shift- 
ing towards high values the range of band fillings where 
order is stabilized. Magnetism, more precisely the pres- 
ence of local moments, strongly reinforces this tendency. 
This should equally apply to other similar alloys, antifer- 
romagnetic coupling becoming competitive at lower band 
filling as in the case of FePt [45] . 

Fruitful discussions with C. Barreteau, J. Los, and G. 
Treglia are gratefully acknowledged. 



of paramagnetic alloys occurs then when is comprised 
between 2.5 and 7.5. However, several alloys, notably 
CoPt with a Nd of 8.3, order rather than phase separate 
for higher values of Nd 0- 

Off-diagonal disorder defined as the difference between 
the band widths of the two constituents Snd = — 
could be relevant here because of its importance when 
mixing Zd and bd elements. Indeed, recently Los et al. 
have shown that the variation of V with N d was ex- 
tremely sensitive to 5nd [43]. Therefore for values of 
5d and 5nd adapted to CoPt, the region of ordered al- 
loys can be shown to be shifted towards higher values 
of N d-, and can include CoPt. Thus, contrary to what 
was argued usually [llj off-diagonal disorder can favor 
ordering instead of phase separation. However, the main 
reason for the ordering of CoPt and several other LIq 
structures seems to be the occurrence of ferromagnetism 
which modifies the relevant number of d electrons [7] . 

The role of magnetism can indeed be illustrated in 
a relatively simple manner in this scheme, by taking 
into account the strong ferromagnetic character of Co-Pt 
(Fig. St and d). The majority spin up band is completely 
full (5 out of the total 8.3 d electrons of the equiatomic 
CoPt are occupied) and does not participate to the co- 
hesion of the system. This results in considering a new 
effective average A^^^ equal to 3.3 d-electrons for a band 
normalized to 5 electrons or 6.6 electrons for a 10 elec- 
tron band. This shift in the effective value of Nd due 
to magnetic considerations pushes the system fully into 
the area oi V > 0, which in turn explains the stability 
of ordered structures and the negative heats of formation 
present in this system. 

Actually it is certainly the presence of magnetic mo- 
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